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SUMMARY 


This paper describes a thoroughly documented experiment that was specif- 
ically designed to test and guide computations of the interaction of an im- 
pinging shock wave with a turbulent boundary layer. Detailed mean flow-field 
and surface data are presented for two shock strengths which resulted in at- 
tached and separated flows, respectively. Numerical computations, employing 
the complete time-averaged Navier-Stokes equations along with algebraic eddy- 
viscosity and turbulent Prandtl number models to describe shear stress and 
heat flux, are used to Illustrate the dependence of the computations on the 
particulars of the turbulence models. Models appropriate for zero-pressure- 
gradient flows predicted the overall features of the flow fields, but were 
deficient in predicting many of the details of the interaction regions. Im- 
provements to the turbulence model parameters were sought through a combina- 
tion of detailed data analysis and computer simulations which tested the sen- 
sitivity of the solutions to model parameter changes. Computer simulations 
using these improvements are presented and discussed. 


INTRODUCTION 


The availability of larger, faster computers, the need to reduce wind- 
tunnel testing, which is time consuming and costly, and the need to provide 
alternate simulation capability for test conditions beyond the reach of prac- 
tical wind-tunnel design have resulted in Increased emphasis on computational 
fluid mechanics. Computations that were not feasible several years ago are 

^Paper to be presented at the Symposium On Flow Separation sponsored by. 
the AGARD Fluid Dynamics Panel at Gottingen, Germany, 27-30 May, 1975, 


377 


now being perforaed routinely. For example: inviscid three-dimensional com- 

putations for speeds ranging from transonic to hypersonic are possible (refs. 

1 and 2); even flows where dominant viscous interactions occur, such as the 
interaction of a shock with a boundary layer, are being computed (refs. 3 
and 4). These interacting flows present the greatest challenge, however, 
because most of the practical applications occur at high Reynolds numbers 
where the flow is turbulent and little is known about turbulence modeling. 

Until recently, most of the techniques for predicting the flow behavior 
in the vicinity of shock boundary-layer interactions were based on experimen- 
tal correlations or approximate solutions to the boundary-layer equations. 

But advances in numerical methods and increased computer speed and capacity 
have resulted in successful attempts to obtain steady-state solutions of the 
complete time-dependent Navler-Stokes equations. Computations of laminar in- 
teractions have been used to illustrate the utility and accuracy of these 
techniques (refs. 3 and 4). Very recent examples for turbulent flows (refs. 
5-7)- Illustrate that such computations using the time-dependent, time-averaged 
Navler-Stokes equations! are also feasible and that they describe the qualita- 
tive features of -the flow interactions. In these examples, closure of the 
conservation equations was accomplished by describing the turbulent shear 
stress and heat flux by eddy-viscosity models and turbulent Prandtl numbers. 
However, a lack of sufficient detailed experimental data precluded verifica- 
tion of these models and efforts to modify them so that the quantitative flow 
features could be more aptly described. . 

The present study was undertaken to provide a formidable first step to- , 
ward understanding the mechanisms that must be modeled before successful nu- . 
merlcal calculations of these complicated flows can be made. It combines ex- 
perimental and numerical methods to guide and verify turbulence modeling for . 
two shock boundary- layer Interaction flows, one. with and one, without separa- 
tion. An axisymmetric experimental arrangement was chosen to. assure purely 
two-dimensional flow (ref., 8). Shocks of two different strengths were im- 
pinged on an established turbulent boundary layer to set up both unseparated 
and separated flows in the interaction zone. Detailed measurements, consist- 
ing of surface-pressure, skin-friction, heat-transfer, and boundary-layer pro- 
files of velocity, static pressure, and temperature were obtained at finely 
spaced Intervals along the surface. Analysis of these data was undertaken to 
define the detailed behavior of the turbulence parameters used to describe the 
shear stresses throughout the interaction regions. Computations were made 
using the time-dependent, time-averaged Navier-Stokes equations employing the 
exact . experimental boundary conditions ‘and algebraic eddy-^vlscosi-ty descrip- 
tions for the turbulent shear stress. The authors reported preliminary prog- 
ress on the separated case in reference 9. 


tlime averaging in these equations is over periods long compared to tur- 
bulence time scales, but short compared. to the time variations of the flow 
field as a whole. 
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NOTATION 


A, A^ Van Driest damping parameter, Eqs. (5) and (6) 
w . 

c specific heat at constant pressure 


S 


specific heat at constant volume, 

local skin friction coefficient, t /(l/2)p u^ 

' ' OO 00 

local Stanton number, q/puc(T-T) 

^ “ p t w 

total specific energy per unit volume, p[c,^T + (y^ + u^) /2] . 


F,G,.H mass-averaged fluxes, Eqs. (1) and (2) 
I intermlttency factor, Eq. (9) 

thermal conductivity • 

mixing length, Eqs. (4) and (20) 
time-averaged pressure 


k 

S, 

P 

Pr 

Pr 


molecular Prandtl number, yc /k 

P 

turbulent Prandtl number 
c 

q^ axial heat flux, -Cp(y/Pr + pe/Pr^)9T/9x 

q^, q , radial heat flux, -c^Cy/Pr + pe/Pr^)9T/9r 
r radial coordinate, distance from model centerline 

t time 

T temperature 

u mass-averaged velocity component in axial direction 

U mass-averaged conservation variables,. Eqs. (1) and (2) 

V mass-averaged velocity component in radial direction 

X axial coordinate, distance from leading edge of shock-wave generator 

y distance normal to model surface 

a shock-wave generator leading-edge angle 
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6 * 


6 * 

c 


£ 

K 

y 

p 

a 


boundary-layer thickness 

boundary- layer thickness at the upstream location of the first measured 
profile station 

kinematic displacement thickness, Eq. (8) 

compressible displacement thickness, Eq. (11) 
eddy viscosity, Eqs. (4) and (20) 

•Von Karman constant, Eq. (5) . ’ 

viscosity 

time-averaged density 


V 9v 3u 
r 3r 9x 


'0 


3v 2 

radial normal stress, p - + -j 

axial normal stress, p - 2y_ + "T y™ 

’ 3x 3 T 

azimuthal normal stress, p - 2y^ ^ "f [r^'*'l^^lx 

total shear stress, Eq. (3) 


V 3v 3u 
r 3r 3x^ 


xr 

Subscripts 

f final axial grid location 

1 Initial axial grid location 

maximum 


max 

o 

r 

t 

T 

X 

w 


location of Incident shock Impingement on cylinder surface In the ab- 
sence of a boundary layer 

radial direction 

stagnation conditions 

total 

axial direction 
wall 

local free stream ahead, of Interaction 
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APPARATUS AND PROCEDURES 


Facility 

The experiment was conducted in the NASA Ames Research Center 3.5-Foot 
Hypersonic Wind Tunnel. This facility operates in a blowdown mode and uti- 
lizes contoured axisymmetric nozzles to achieve a uniform Mach number, and an 
open-jet type test core approximately 0.7 m in diameter and 4 m in length. 

The present tests were all performed at a nominal free-stream Mach number of 
7.2 and with nominal values of total temperature and pressure of 695 K .and 
34 atmospheres, respectively. The corresponding nominal value of free-stream 
unit Reynolds number was 10.9 x 10^ m“^ and the useful test time was about 3 
min. 


Model 

A cone-ogive cylinder, 330 cm in length and 20.3 cm in diameter was used 
as the test surface (see fig. 1). An annular shockwave generator, 51 cm in 
diameter, mounted concentric with the cylinder axis was used to generate shock 
waves of two different strengths by beveling the sharp leading edge at either 
7.5° or 15°. The generator could be translated in a direction parallel to the 
cylinder axis so that the entire Interaction region could be passed over se- 
lected survey stations. 

Interchangeable instrumentation ports, 12 cm in diameter and specifically 
contoured to match the cylindrical surface, were located at 25 cm Intervals 
along the cylinder in a single line, and every 50 cm in another line 180° 
around the body. One port was Instrumented with a floating element skin- 
friction balance. Another was used to accommodate either pitot and static 
pressure probes or total temperature probes. The probes were positioned by a 
mechanism contained inside the cylinder and automatically actuated from out- 
side the tunnel test section. The remaining ports were Instrumented with 
thermocouples spot-welded to the inner surface every 1.25 cm and with static 
pressure taps. Static pressure taps were also located every 5 cm along the 
entire cylinder between ports. 


Test Procedure 

Data were collected from a series of tests with the tunnel operating at 
the nominal conditions described above. In separate tests without the gener- 
ator it was determined that a fully developed turbulent boundary layer with 
negligible axial pressure gradient was established over the cylinder surface 
between 100 and 300 cm from' the model tip (ref. 10). With the generator in 
place, the nominal measured boundary-layer parameters ahead of the interaction 
for the 7.5° and 15° generator angles, respectively, were: edge Mach number, 

6.7 and 6.9; boundary- layer thickness, 3.2 cm and 2.7 cm: and Reynolds number 
based on boundary-layer thickness, 0 . 23 x 10® and 0.2 x 10®. The model wall 
temperature was essentially constant at a value of 300 K. 
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Surface pressure, skin friction, and heat transfer were obtained at small 
intervals by moving the shock-wave generator in an axial direction durli^ the 
tests. The boundary-layer thickness Increased about 10% over the distance of 
25 cm which corresponded to the difference* between the farthest upstream and 
downstream positions of the generator. The difference in boundary-layer 
thickness had little influence on the results, provided they were compared at 
equivalent distances from the leading edge of the generator. 

* Skin friction was measured with a contoured floating-element balance 
whose sensible element was 0.95 cm in diameter. Calibrations of the gage be- 
fore and after each test run were repeatable to within 5%. The skin-friction 
data were corrected for buoyancy effects resulting from the axial pressure 
gradient. Corrections were less than 10% of the measured values, except in 
the regions of minimum skin friction for the 7.5“ generator tests and near 
separation and reattachment for the 15“ generator tests, where they were as 
high as 50% of the measured upstream zero pressure gradient values. The heat- 
transfer rate was measured using the thln-wall transient technique. Longitu- 
dinal conduction errors were computed and found to be less than 5% of the 
measured rates and so no corrections were applied to these data. 

Velocity, density, and pressure profiles were obtained from pitot and 
static pressure and total temperature surveys. Each survey was taken during 
a single test run, and its axial location was established prior to the run by 
prepositioning the shock generator. In the interaction region, surveys were 
obtained every 2 cm for the 7.5“ generator tests and every 2.5 cm for the 15“ 
generator tests. Downstream of the interaction the corresponding distances 
between survey stations were Increased to 4 cm and 5 cm, respectively, for the 
two generator angles. In the reversed flow region established with the 15“ 
generator, pitot measurements were obtained in upstream and downstream direc- 
tions to help establish the extent of separation. When traversing the bound- 
ary layer, the probes were stopped at each location for a few seconds to avoid 
time, lags in the measurements, and static pressures at the model surface were 
monitored continuously to verify interference-free data. 

To verify that the model was aligned with the free-stream flow direction, 
surface-pressure measurements at selected axial positions were obtained at 90° 
intervals around the model, and skin-friction measurements at selected axial 
positions 180“ apart. Comparisons of these data around the model showed vari- 
ations that were wlthiln the experimental accuracy of the measurements. For 
the 15“ generator tests, separation and reattachment lines around the model 
were also measured using a surface oil film technique. The results, verified 
an axlsymmetrlc separation zone. 

A more complete description of the test procedure and data accuracy along 
ffith tabulations of all the test data are given in reference 11. 

GOVERNING EQUATIONS AND NUMERICAL PROCEDURES 

The equations and numerical procedures were first presented by the au- 
thors in reference 9. For completeness, some of that information is presented 
again in this section. 


382 



Equations and Boundary Conditions 

The mass-averaged Navierr-Stokes ^equations for compressible flow,^ex~- .. 
pressed in cylindrical coordinates with axial symmetry assumed, were used to 
predict the flow throughout the interaction region. The equations and. the. . 
concept of mass averaging are discussed in reference 12. The turbulent Rey-. 
nolds stress and heat-flux terms in these .equations are related to. the mean 
flow gradients of velocity and temperature by eddy-transport coefficients that 
are added to the molecular-transport . coefficients . Additional, restrictions on 
the equation system include the perfect gas assumption, constant specific . . , 
heats, the Sutherland viscosity iaw,. and zero bulk viscosity., The resulting , 
equations are ... . , , . . .. 


^ 

3t • 3x • 3r 


p 


( ' ' \ 
pu 

pu 

F = r 

pu^ + a 

X 

pv 


puv + T 
; xr 

e 


(e + o )u + T v3+ q 



X , xr . ■X.-' 

'O' 


' pv 

0 


puv ■+ T 

xr 

^6 

' G * r. 

pv^ ^'r 

0 


(e + 0 )v + T u + q 

V > 





( 2 ) 


Figure 2 shows the computational domain. The conditions on the upstream 
boundary, were prescribed by a combination of an inviscid, method of character- 
istics program (ref. 13), and, a boundary-layer program (ref. 14) modified for 
turbulent flows by Marvin and Sheaf fer. At the upstream boundary position, 
ithe experimental arid computed ihciderit shock waves were aligned and the ' 
boundary- layer program was run for an x distance that ensured a match of 
experimentally and numerically determined displacement thickriessesl Along the 
cylinder surface, r = r^, or y = 0, the boundary conditipns used were the 
viscous, rio-slip conditions, u = v = 3p/3r =0, T = T^, while along the outer 
boundary they were the iriviscid, free-slip conditions, v = 3p/3r = 3u/3r 
= 3T/3r =0. At the downstream boundary the derivatives of all variables 
were set to zero, e.g., 3u/9x = 0. The initial conditions within the com- ■ 

putatiorial domain used to start the solutions were obtained by setting^ the 
values of all variables equal to theiir inflow boundary values at the same 
vertical station, that is, f(x, r, 0) = f(x^, r, 0), x^^ < x ^ Xf. To 
restart solutions or make modifications to turbulence parameters, initial val- 
ues of the variables were set equal to their computed values obtained during 
the last time step of the previous solution.' . < =. , . 
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Numerical Procedures 


The finite difference scheme used to solve equation (1) is the same as 
that developed originally by MacCormack (ref. 3) and applied more recently to 
two-dimensional turbulent shock boundary- layer .interactions by Baldwin and 
MacCormack (ref. 15). The numerical techniques employed in the present study.,, 
along with all the pertinent difference equations, and some special proce- 
dures are reported in reference 9. The computational domain was subdivided 
into four .subgrids with each subgrid divided into a number of uniform grid 
cells with spacing Ay. Finer spacing was employed near the wall. A total of 
78 cells . in the y direction, was used; uniformly spaced grid with 48 cells in 
the . X direction was used with spacing ,.Ax = 0.80 or 0.635 cm for the 7.5° and 
15° cases, respectively. Some of the solutions presented later coyer axial 
distances greater than those obtained with the 48 cells in the x direction. 
Those solutions were achieved by redefining the upstream boundary to coincide 
with a position about 3 cm ahead of the downstream boundary from converged 
solutions and then continuing the solutions on downstream for another 48 
cell points in the x direction. 

The solutions were advanced in time following the procedure described in, 
reference 9. Steady-state convergence was assumed when solutions from at 
least 20 successive tlm^ steps showed little or no change. Computation times 
to achieve these fully converged solutions on a CDC 7600 were 3-4 hrs for the 
15° generator cases and about 1 hr for the 7.5° cases. 


RESULTS AND. DISCUSSION 

. . •'’j. 


Experimentally Determined Flow-Field Features 

Figure 2 depicts the major features .of the shock-wave, ’ boundary-layer 
Interaction zone. The sketch is based on boundary-layer survey measurements 
and shadowgraphs taken during the experiments using both the 7.5° and 15° 
shock-wave generators. The incident shock wave, weakened and curved somewhat 
by the expansion fan emanating from the corner formed by the leading edge and 
the body of the shock-wave generator, impinges on the incoming boundary layer. 
The subsequent increases in surface pressure cause the boundary layer to 
thicken, or even separate in the case of the strongest Incident shock wave, 
and induce a shock wave. Thereafter rapid flow turning and boundary-layer 
thinning occur and a recompressiqn shock is formed. 

Figure 3 shows the surface measurements obtained for the two shock-wave 
generator angles. Surface-pressure, skin-friction, and heat-transfer coeffi- 
cients are shown as functions of a normalized interaction distance centered 
about Xq, the location of the Intersection of the inviscid incident shock 
wave with the body surface in the absence of a boundary layer. With the 7.5° 
generator, the pressure rises continually through the interaction; the skin 
friction d.ecreases initially in the presence of the adverse pressure gradient 
and rises thereafter in the recompression region where the boundary layer 
thins. The heat transfer follows a behavior similar to that of the pressure. 
No separation was observed in this case, either from the skin-friction mea- 
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surements or from oil-flow patterns that were established and photographed 
during some special tests. In contrast, with the 15° generator the initial 
Increase in pressure levels off in a plateau, and negative values of skin 
friction were measured, both characteristics typically associated with separa- 
tion. Rapid Increases in pressure and skin friction occur downstream of re- 
attachment. The heat transfer rises continually through the interaction until 
the surface pressure decreases. For both generator angles, the decay in pres- 
sure, skin friction, and heat transfer downstream of the Interaction is a di- 
rect result of the expansion fan emanating from the corner of the generator. 

Figures 4 and 5 present constant static pressure and velocity contours 
and Illustrate further the details of the two interaction cases under investi- 
gation. The contours were constructed from the profiles of velocity and stat- 
ic pressure across the boundary layer obtained from pitot and static pressure 
and total temperature measurements taken at small Ax and Ay intervals. Com- 
plete tabulations of these profile data can be found in reference 11. Loca- 
tions of the incident, induced, and recompression shocks are easily recognized 
in the pressure contours. Note that for the 15° generator the pressures are 
higher than the wall-peak pressure in the compression region downstream of the 
intersection of the incident and induced shocks. The velocity contours illus- 
trate more dramatically the differences between the two “"flow cases. For the 
larger generator angle the flow velocity near the wall in the vicinity of the 
interaction is highly retarded and achieves negative values associated' with 
the reversed flow in the separated zone. 

A precise determination of separation and reattachment points for the 
separated case was made difficult because of the unsteady nature of the separ- 
ated flow and the relatively large diameter of the skin friction element. A 
detailed discussion of this unsteady phenomenon is given in reference 9. A 
best estimate of the extent of the time-averaged separation region was obtain- 
ed from data obtained with forward and backward facing pitot tubes. These 
data^ obtained at fixed values of y with x varied by moving the shock gener- 
ator, indicated a separated region extending from (x-Xq)/6q = -3.15 to -1.68. 
This region is somewhat larger than the skin-friction measurements indicated. 
Locations of the separation and reattachement points from the pitot measure- 
ments are shown on the abscissa of the skin-friction plot. 

Numerically Simulated Flow-Field Features 

As previously mentioned, the turbulent Reynolds stress and heat-flux 
terms in equation (2) were assumed to be related to the mean flow gradients of 
velocity and temperature by algebraic eddy-transport coefficients that were 
simply added to the molecular- transport coefficients. Mainly, this choice was 
dictated by considerations of economy in the computer program. Although such 
a model may be restrictive in its application to other new flow situations, it 
suits our current objectives of (1) defining the eddy viscosity field that re- 
sults when a shock impinges on a turbulent boundary layer, and (2) determining 
whether improvements in the eddy viscosity description used in the numerical 
simulations can be made by a close examination of the experimental data. 
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Results using a baseline turbulence model = — The first numerical simula- 
tions were obtained with a modified, two-layer, Cebecl-Smlth (ref. 16) eddy 
viscosity model. Modifications suggested byr Cebeci to account for pressure 
gradient effects were not used. The shear was expressed as 


T = = (y + pe) 
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In the inner layer the following mixing length description for the eddy vis- 
cosity was used: 
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with the Von Karman constant k = 0.4 and the Van Driest constant A =26. 
In the outer layer, the eddy viscosity was given, following Clauser, oy 
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The above definitions of the displacement thicknesses differ somewhat from the 
conventional ones. The present choice was dictated by the fact that, in the 
first stages of developing the Navier-Stokes code for the shock interaction 
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problem, overshoots in the velocity profiles during early time steps made it 
difficult to select the edge of the viscous layer. Therefore, the value of 
ymax simply taken as the boundary- layer thickness ahead ' of the interac- 
tion, and (pu) mav were taken as the local maximum values between 

the wall and Ymax* boundary between the two layers was determined by the 

value of y where 

■' inner outer 


The heat flux was expressed in terms of the eddy viscosity by 



( 12 ) 


where Pj.^ = 0.9. 

Results of the computations using this baseline turbulence model are com- 
pared with the experimental data in figures 6, 7, and 8. Overall, the compu- 
tations predict the qualitative features of the two flows remarkably well con- 
sidering the simplicity of the turbulence model, but a closer examination of 
the comparisons points out the major limitations of the computations. 

The overall surface pressure .rise for the flow with the 7.5° shock-wave 
generator (fig. 6a) is predicted reasonably well, except for the location of 
the initial rise in pressure. The corresponding predicted rises in both skin 
friction and heat transfer lag the. data in the interaction region, reflecting 
the Inability of the s:mple turbulence model to predict any upstream influ- 
ence, but the final predicted levels downstream agree reasonably well with the 
data. The skin friction prediction shows separation at the surface but the 
measurements do not. The pressure contour comparisons (fig. 7a) show that 
the computation predicts only the incident and reflected shocks whereas the 
experimental data show the presence of an induced shock. The comparison of 
the streamline contours (fig. 8a) shows the prediction of a zero velocity 
line just off the surface accompanied by a small region of reversed flow which 
is not present in the experimental data. 

4 

Similar conclusions can be made from comparisons of the numerical compu- 
tations and the data for the flow with the 15° shock-wave generator (figs. 6b, 
7b, and 8b). With the baseline turbulence model, no upstream influence is 
predicted and in this separated flow case no plateau pressure is predicted. 

The induced shock wave caused by the large separation in the experiment is not 
predicted because the computation predicts such a small separation height. 

For this separated flow case, the baseline turbulence model results in good 
prediction of the overall pressure and skin friction rises and their subse- 
quent decay; but the heat transfer is substantially under predicted (see fig. 
6b). 


Baseline model modifications . — Attempts were made to guide changes in the 
turbulence model by combining data analysis and trial and error solutions in 
the actual Navier-Stokes code. Ideally, these changes could have been guided 
entirely from data if absolutely reliable shear-stress measurements had been 
available. But, since attempts to directly measure the shear’ stress through 
these interactions have so far produced unsatisfactory results, the boundary- 
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layer profile data were used. A significant degree of uncertainty is intro- 
duced in this procedure, however, because the inertial forces dominate the 
momentum balance, especially in the outer portions of the flow. Conversely, 
relying solely on trial and error solutions to the code itself to guide model- 
ing changes would be time consuming and perhaps unsuccessful if attention was 
not confined to physically meaningful changes. 


The shear-stress and heat-flux distributions through the boundary layer 
were evaluated by the use of experimental profile data to solve the following 
equations based on the boundary- layer approximation: 
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By performing integration with respect to y before differentiation with re- 
spect to X and by employing the conservative form of the variables, for ex- 
ample (p + pu^) , it was expected that errors in the momentum and energy bal- 
ances could be minimized. Despite these precautions, not all ..the shear pro- 
files approached zero at large distances from the wall where they should have. 
In some of these cases it was possible to adjust the inertial balance across 
the boundary layer so that zero shear was achieved at the edge of the thermal 
boundary layer. These adjustments were usually small for the 7.5° shock wave 
generator profiles, but somewhat larger for the 15° generator profiles. 


Figure 9 shows the shear profiles resulting from these momentum and ener- 
gy balances for axial locations ahead, within, and downstream of the two in- 
teraction regions. For the flow with the 7.5° generator the maximum shear 
stress within the boundary layer builds up rapidly within the interaction re- 
gion as the adverse pressure gradient increases; after peak pressure is 
reached this maximum shear relaxes toward its initial level but at a very slow 
rate. At the farthest downstream location where the pressure gradient is 
favorable the shear near the wall decreases and subsequently Increases, indi- 
cating that the shear in the outer extremes of the boundary layer is adjusting 
to the local flow gradients more slowly than the shear near the wall. Similar 
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conclusions can be reached for the separated case using the 15° generator. In 
this case, data were available for a larger downstream interaction distance, 
and the shear appears to have adjusted to the local flow gradients. Also, for 
this separated case, the maximum shear in the boundary layer continues to in- 
crease downstream of reattachment and slightly beyond the location of peak 
pressure. 


A maximum mixing length was determined at each of the profile survey sta- 
tions by dividing the experimentally deduced shear distributions by the mea- 
sured velocity gradients, plotting the results, and choosing the maximum value 
of mixing length. For those cases where the mixing length continuously in- 
creased with distance from the wall, the value of maximum mixing length was 
chosen at the point where the first significant departure from a linear mixing 
length distribution occurred. These maximum values are shown in figure 10 as 
a function of the interaction length parameter. The extremes on the bars rep- 
resent the uncertainty introduced by using shear profiles evaluated either 
directly from the momentum balances or from momentum balances modified by ad- 
justing the inertial terms to Insure zero shear at the edge of the thermal 
boundary layer. The uncertainty was largest in the region downstream of reat-' 
tachment for the 15° generator case. In both cases the maximum mixing lengths 
tend to decrease in regions of adverse pressure gradient, where the boundary 
layer thickens, and to increase in regions of favorable pressure gradient. 

Near the surface where the Importance of the inertia and convection terms 
in the momentum and energy balances diminished, attempts were made to evaluate 
K and A^. The shear profiles at each survey station were analyzed by inte- . 
grating the following system of equations to obtain values of velocity and 
temperature as a function of y out to distances where the estimated errors 
in the Inertial balance became significant compared to the magnitude of the 
local shear (usually this consisted of about 9 measured points away from the 
wall) : 
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Best fits to the velocity and temperature profile data near the wall were 
achieved by repeated integration of these equations until thfe sum of the root 
mean square of the differences between the predicted and measured velocities 
and temperatures was minimized. The procedure was automated for solution on 
a CDC 7600 and initiated by inputting shear and temperature profiles from the 
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momentum and energy balances along with the measured values of wall shear. An 
optimization routine (described in ref. -16) was initiated and values of k; and 

and were sought to satisfy the minimization criteria, (Tliis optimiza- 

tion routine was developed by Garret N. Vanderplaats , of NASA Ames, who helped 
to ' implement it for the present application.) The root mean square of the 
residuals never exceeded 5% and in most cases was less. Exceptions to this 
residual band occurred for the 15° generator case in the separated) region., ;> - 

However, it was still possible to achieve this band at these profile stations, 
provided the input wall shear was also, considered part of the optimization 
routine along with the other parameters. This was not surprising, however, 
since accurate skin-friction measur^ents were difficult to make in this sep- 
arated region. Examples of the best fits to the velocity and temperature 
data hear the upstream edge of the separated region are shown as the solid 
lines in figure 11. Two important aspects of this example are noteworthy. 
First, at this station, the wall shear needed to achieve a best fit was 
-5.28 N/m2 whereas the direct measurement was 12.1 N/m2. Consideringrthafe v; 
the station is near the separation point where experimental accuracy is poor 
and' 'that the separation point is unsteady i this disparity is not unreasonable. 


Second, at this station, values of 


are significantly lower than 


their corresponding undisturbed values, .0.4 and 26- ‘ • i. i. if 

Figure 12 presents the values of k arid A^ required to achieve these’ 
best fit velocity and temperature profiles for both of the .interaction cases 
being' studied. The error bands on the symbols again represent the uncertainty '- 
Introduced by using shear profiles evaluated either directly from -the momentum’ * 
balances or from momentum balances modified by adjusting the inertial terms to 
Insure zero shear at the edge of the thermal boundary layer. Thd results for 
the 7.5® generator show that A^ decreases" in the vicihit;y of the interaction'* 
while K increases. ' These results We interpreted as indicating a hiW^i^ 
eddy viscosity throughout the region' than that predicted by the baseline tur- 
bulence model. This could explain why the numerical simulations using the 
baseline turbulence model predicted separation. For the 15° generator case 
where separation was present; similar trends in k and A^ are apparent. How- 
ever, the uncertainties within the separated region precluded any precise de- 
termination of the parameters. Appaxehtly," both k arid A^ decrease ahead of 
and in the upstream portions of the sepafatW region. At the downstream edge 
of the separated region near reattachment; ic has increased considerably Wd 
A^ is also increasing. The physical iriterpretatlon of these, results suggests ‘ 
that .ahead of separation the sublayer of" "Wner region has a somewhat higher 
eddy viscWity' than the' Wseline moWl would pfWict," but' in the logarithmic ‘ 
region the viscosity is somewhat lower than the baseline model would predict. 

In the separated region, an interpretation is more difficult to arrive at, but 
apparently, at least near the reattactoent point, the eddy viscosity is higher 
across the sublayer,’ and logarithmic regions than would be predicted by the 
baseline model. i 

The turbulence model mixing length formulation suggested by the foregoing 
data analysis was introduced into the Navier-Stokes computer code in the fol- 
lowli^ way; ’ 



In the ''inner region • . . i ? > 



(19) 

( 20 ) 

. .WV : / 

' K '{ ^ ''i I 

( 21 ) 





Employing ^^niax this model formulation eliminated the need for . arbitrarily 
defining compressible and incompressible displacement .thickness as was. the 
case with. the baseline model. The boundary between the two regions is' deter-, 
mln^ by.,the value of y where i =' The heat flux equation remained 

the , same, as for the baseline model (see eq. ,12). . . ! 

• • •• , • .* . f. . ' • . • , ‘ 

.**}, .> • ‘ • . . ; • ‘ ’ • • • * . 

.. Computer simulations were next obtained using the experimental data anal- ' 

ysls as a guide for evaluating the parameters it ina v(x) , ic(x) , and. ^(x) . . 

First, values directly from the data analysis were used. Examination of the 
resulting computer simulations made it apparent that adjustments to the param- 
eters would be needed before the simulations would predict the experimentally 
determined features of the flow fields in the interaction region. Therefore, 
a trial and error procedure was initiated to arrive at more appropriate dis- 
tributions of the turbulence parameters. , The procedure is still underway at 
this. time. Before discussing the results of this procedure to date, some in- 
teresting observations can be pointed out. The simulations were all more ' 
sensitive .to modifications of the inner region model parameters than the qutjer .' 
maximum mixing length parameter. In tHe inner region itself, A^(x) modifica- 
tions tended to affect the solutions' more than those for k(x), especially in 
the separated case. However, best results have, been achiey,ed .with modif lea- ,\“ 
tlon to both A^(x) and k(x) . 

Results using a modified turbulence model . — The turbulence parameter var- 
iations used in the latest computer simulations employing the complete Navier- 
Stokes equations are shown as the solid and dashed lines in figures 10 and .12. 
The maximum mixing length variation (fig. 10), employed in all simulations, 
corresponded to the mean variation exhibited by the data analysis. In the in-r 
ner layer region (fig. 12), values of and k had to be altered to obtain 
better predictions of the experimental surface and flow-field data. For the 
7,5° generator case these variations in and k follow the trends exhibited 
by the data analysis. The largest disparity between the values used in the 
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simulations and those deduced from data "analysis is in • at the beginning 
of the interaction region. As noted before, decreasing A^ and increasing 
both result in a larger eddy vlscositysin the inner layer region. For the 15 
generator case the variations in the '‘inner layer turbulence parameters- follow 
the same trends as the data analysis except near reattachment ; Part of the 
reason for this disparity is that the- data. analysis in this region was com- 
pleted ’only recently and , there was insufficient . time to. .obtain,, converged ^simu- 
lations before preparing this manuscript. rij> - try 


The results of the computer simulations using the turbulence parameter 
variations described above are presented in figures 6, 7, and 8. For the 7.5*’ 
generator case'there' is •obvious; improvement, over,' -the; baseline .predictions.- of 
■’ “sWin friction'and' heat transfer inithe.tlntefaetion.'regionr,i(see:,fjLg.i'6a)> * . Sep- 
"aratioh’ is not predicted and 'the" upstream, influence, of the interact ipnocoln-: 

'■ “cides with- that^bbserved ’from" the data'r^ ^Downstream of.. the peak pressurej.lpca- 
' tion the heat- transfer predlctlbh-'ls ;not ; as rgood as the -baseline ’modeljpr.edic- 
tions," but this could be explained -by"' an incorrect cho.ice .of a'.constant' turbu- 
i'ent ‘Prandt-l number ; ■ At ’ this stage of-’'model development, this disparity ■ ini- 
heating prediction is hot considered^ crucial.* because the solutions of the mp- 
menitum and ' energy equations are loosely^ coupled ., The turbulence model . changes 
had little effect on the surface pressure prediction. .The pressure and veloc- 
ity contours (figs. 7 and 8) are not changed significantly from those for the 
baseline model, except that no reversed flow region is predicted with the mod- 
ified model. The reflected shock observed in the data is still not predicted. 
The main reason for this is that the numerical simulations fail to show a sig- 
nificant thickening of the boundary layer at the start of the interaction re- 
gion. ‘ " - i'' 

Results using the turbulence model modifications for the 15 generator 
case are also shown in figures .6, 7, and 8. The predicted separation bubble 
size* increased considerably with a .corresponding prediction of upstream influ- 
ence and a plateau in the surface pressure. The pressure contours show the 
presence of an induced shock wave similar to that observed in the experiment. 
Obviously-, substantial deficiencies still exist in the predictipn using the 
modified model. The separated bubblensize; is still smaller in height than the 
experiment indicates; the plateau pressure; is only about half the measured 
value; and all three surface quantities show substantial differences with the 
measurements in the interaction region. It is felt, however, that a signifi- 
cant improvement' in the surface predictions can be made by including the lat- 
est data analysis values of A^ and k near reattachment. This should 
shorten the extent of predicted separation and shift the rise in skin-friction 
and heat-transfer upstream. * : - • 


• ' CONCLUDING. REMARKS 

A detailed experimental investigation of the mean flow throughout two 
shock-wave boundary- layer interaction' regions, one with separation and one 
without, has been presented. Although the interactions were very complex, the 
mean data were of sufficient detail and quality to assess the validity of nu- 
merical simulations and to guide turbulence model changes. 


O A 


• Numerical solutions, employing the full time-averaged Navier-Stokes 
equations along with algebraic eddy-viscosity models appropriate for zero- 
pfressur e-grad lent flows, predicted the overall features of the flow fields, 
blit they were seriously deficient in predicting the details of the interaction 
regions. Through a combination of data analysis and trial and error computer 
simulations, which tested the sensitivity of the solutions to turbulence model 
parameter changes, the agreement between numerical predictions and experiment 
was improved. . 

Although the improvements fell short of identifying an optimum model, 
several important trends regarding the two-layer, algebraic eddy-viscosity 
model- can be noted. The inner layer model parameters had substantially more 
Tliifluence- on -the numerical simulations than the outer layer parameter. Where 
the 'boundary layer was unseparated, the teddy viscosity in the inner layer re- 
gion-had -to be increased substantially over that predicted by a zero-pressure- 
giradient, two- layer model; otherwise separation was predicted. Therefore 'the 
•simple' zero-pressure-gradient model cannot be used to predict locations, of 
separation. Where separation did occur, the results. were less clear, but the 
eddy viscosity in the sublayer region had to be increased while in the logar- 
'ithmlc region it had to be reduced somewhat; otherwise the separation bubble 
'Size was substantially underpredicted.. 
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Figure 3.- Measurements along the model surface 
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Figure 4.- Static pressure contours Figure 5.- Velocity contours obtained 
obtained from flow-field measurements. from flow-field measurements. 
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Figure 6.- Comparisons of computations and measurements along the model surface 
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Figure 9.- Shear profiles for various 
axial locations along the model. 
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Figure 10.- Axial variation of 
maximum mixing lengths. 
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Figure 11.- Comparison of experimental 
and best fit velocity and tmperature 
profiles employing a modified Inner 
eddy-vlscoslty model. 
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Figure 12.- Axial variation of 
A^ and k. 




